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SUMMARY 


We  have  examined  the  theory  of  both  electrical  and  thermal 
breakdown  in  explosive  materials.  Application  of  thermal  breakdown 
theory  had  to  be  modified  to  account  for  the  reaction  kinetics 
while  electrical  breakdown,  whose  onset  occurs  at  the  steady 
state  temperature,  was  found  to  follow  breakdown  theory  in  inerts. 

Use  of  breakdown  theory  to  predict  critical  fields  in 
explosives  seems  feasible  but  is  hampered  by  the  lack  of  an 
adequate  data  base.  This  lack  of  experimental  data  also  prevents 
a  resolution  of  the  question  of  possible  electric  field  effects 
on  thermal  explosion  parameters.  It  should,  however,  be 
relatively  simple  to  distinguish  between  thermal  and  electrical 
breakdown  in  any  given  experiment  by  a  consideration  of  initial 
temperature  and  explosion  induction  time. 


ELIHU  WMET 
By  direction 


Accession  For _ 

II  I  I  S  GRA&I 
DIIC  TAB 

Unannounced 
Justification — 


iii  * 

list  r  5  V’t'- — I 
Avail?.*-- *  !.V«  _ 


.  ’  o  r 

o-c—1 


If! 


NSWC  TR  79-331 


CONTENTS 

Page 


I.  INTRODUCTION . 4 

II.  ELECTRICAL  BREAKDOWN . 4 

1.  Formulation  of  Breakdown  Criterion . 4 

2.  Electron-Phonon  Interaction . 7 

a.  Polar  Molecules. . 7 

b.  Non-Polar  Molecules . 9 

3.  Relaxation  Times.  . . 11 

4.  Rate  of  Energy  Transfer  to  Lattice . 14 

5.  Calculation  of  Critical  Fields . 15 

6.  Time  Lag  Before  Breakdown . 15 

7.  Characteristic  Results . 22 

8.  Critical  Electron  Density . 22 

9.  Effects  of  the  Electron  Distribution . 24 

10.  Application  to  Explosives . 26 

III.  THERMAL  BREAKDOWN . 26 

1.  Thermal  Breakdown  in  Inerts . 26 

2.  Thermal  Breakdown  in  Reactive  Media . 27 

3.  Crank-Nicholson  Method . 2G 

4.  Numerical  Integration  Results  and  Discussion.  .  .32 

5.  Influence  of  the  Conductivity  Work  Function.  .  .  32 

6.  Thermal  Breakdown  Characteristics . 36 


IV.  SUMMARY  AND  CONCLUSIONS 

V.  ACKNOWLEDGEMENT.  .  .  . 

BIBLIOGRAPHY . 38 

APPENDIX  A  CALCULATION  OF  THE  INTERACTION  CONSTANT.  .  .  .  A-l 

APPENDIX  B  COMPUTER  PROGRAM  FOR  SOLUTION  OF  THE  THERMAL 

BALANCE  EQUATION  UTILIZING  A  CRANK-NICHOLSON 
SCHEME . B-l 

APPENDIX  C  THE  TRIDIAGONAL  MATRIX  AND  THOMAS  ALGORITHM.  .  C-l 


2 


NSWC  TR  79-331 


ILLUSTRATIONS 


Figure 

1 

2 

3 

4 

5 

6 

7 

8 


Critical  Breakdown  Field  as  Function  of 
Temperature  (°K)  in  Polar  Crystals, 

Critical  Breakdown  Field  as  Function  of 
Temperature  (°K)  in  Non-Polar  Crystals 
Graphical  Solution  of  Equation  (68)  for 
the  Collision  Generation  Number  n 
Grid  for  Finite  Difference  Computation 
m  a  Crank-Nicholson  Scheme  . 

Computational  Geometry  -  A  Semi-infinite ’slab 
temperature  Evolution  for  Different  Bath 
Temperatures  TB  . 

Temperature  Evolution  for  Different  Bath* 

Temperature  Tg . 

Dependence  of  Critical  Temperature  for 
Thermal  Explosion  on  the  Conductivity 
Work  Function .  * 


Page 


16 

17 

21 

29 

30 

33 

34 

35 


3 


NSWC  TR  79-331 


I .  INTRODUCTION 

Dielectric  breakdown  is  a  universal  response  of  materials 
to  the  application  of  sufficiently  intense  electromagnetic 
fields.  While  the  interaction  of  electromagnetic  radiation 
with  ine^t  dielectrics  has  received  attention  from  both 
theorist^"6  and  experimentalist7-9,  relatively  little  work 
has  been  done  on  the  coupling  of  radiation  to  explosive  materials, 
a  matter  of  some  recent  interst.  We  may  also  note  experimental 
evidence  indicating  generation  of  intense  polarization  fields 
in  a  detonating  explosively-^2.  An  understanding  of  breakdown 
phenomena  may  thus  help  elucidate  fundamental  detonation  processes, 
even  in  the  absence  of  externally  coupled  fields. 

Theories  of  breakdown  in  solids  divide  broadly  into  studies 
of  thermal  breakdown  and  purely  electrical  breakdown,  with 
a  number  of  subcategorj.es  for  the  latter.  We  shall  briefly 
review  each  process  and  consider  the  application  to  an  explosive 
medium. 


II.  ELECTRICAL  BREAKDOWN 

1.  FORMULATION  OF  BREAKDOWN  CRITERION.  In  electrical 

breakdown  we  assume  that  the . electron  distribution  depends 

directly  on  the  applied  electric  field  t.  The  critical  breakdown 

iFrohlich,  h7,  "Theory  of  Electrical  Breakdown  in  Ionic  Crystals," 
Proc.  Roy.  Soc.  A160,  230  (1937). 

2Frohlich,  H.  and  Paranjape,  V. ,  "Dielectric  Breakdown  in  Solids," 
Proc.  Phys.  Soc.  Lond.  B69,  866  (1956). 

20'Dwyer,  J.  J. ,  "Dielectric  Breakdown  in  Solids,"  Adv.  in  Phys. 

7,  349  (1958). 

■Seitz,  F. ,  "On  the  Theory  of  Election  Multiplication  in  Crystals," 
Phys.  Rev.  76,  1376  (1949). 

^Stratton,  J.  A.,  "Theory  of  Dielectric  Breakdown  in  Solids," 

Progress  in  Dielectrics,  Vol.  3,  Ed.  Birks  (1961,  Wiley). 

6Von  Hippel,  A.,  "The  Electrical  Breakdown  Strengths  of  Ionic 

-Crystals  as  Functions  of  Temperature,"  Phys.  Rev.  56,  941  (1939). 

'Vorobev,  A.,  "Anisotropy  of  the  Dielectric  Strength  of  Rocksalt," 
Soviet  Phys.  Solid  State  4,  1441  (1962). 

®Caspari,  M. ,  "Direction  of  Breakdown  in  Alkali  Halide  Crystals," 
Phys.  Rev.  98,  1679  (1955). 

9Davisson,  J.  W. ,  "Directional  Breakdown  in  Crystals,"  Progress 
in  Dielectrics  Vol.  1,  Ed. -Birks,  (1959-Wiely) . 

lOHayes,  B 7 ,  "iTte  Detonation  Electric  Effect,"  J.  App.  Phys. 

38,  507  (1967). 

HDremin,  A.  et.al,  "Shock  Induced  Electrical  Polarization," 

Proc.  6th  Int.  Symp.  on  Detonation,  (1975) . 

■^Yakushev,  v. ,  et.al,  "On  the  Measurement  of  Polarization  Relaxation 
Time  in  Homogeneous  Explosives,"  Zh.  Experm.  Teor.  Fiz.  54, 

396  (1961). 
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field  Fq  will  be  chat  marking  the  onset  of  some  instability 
in  the  conduction  current. 

We  consider  first  the  case  of  low  carrier  density  and 
follow  the  average  behavior  of  a  single  electron  interacting 
with  the  applied  electric  field  and  with  the  lattice.  As  the 
density  increases,  electron-electron  interactions  become  important 
and  a  single  electron  picture  will  no  longer  suffice.  We  will 
later  calculate  the  critical  density  for  a  single  electron 
approximation  and  consider  the  necessary  modifications  when 
the  low  density  restriction  is  violated. 

For  a  steady  state  to  exist,  an  equilibrium  must  be  maintained 
between  the  energy  supplied  to  the  electron  by  the  external  field 
A(F,E) ,  and  the  net  energy  transferred  from  the  electron  to  the 
lattice  via  collisions  with  phonons,  which  we  shall  designate 
B(E).  E  represents  the  total  electronic  energy. 

We  consider  an  applied  electric  field  F  pointing  along 
the  z-axis  of  our  cartesian  coordinate  system.  A  conduction 
band  electron  will  be  accelerated  by  the  field  feeling  a  force 


dpz 

3t~ 


FIELD 


eF 


(1) 


We  adopt  the  convention  that  a  previously  defined  vector  symbol 
appearing  without  an  arrow  refers  to  the  absolute  value.  The 
buildup  of  momentum  in  the  Z-direction  will  be  retarded  through 
collision  with  the  lattice  leading  to  an  average  drag  force 


dpz 

dt 


LATTICE 


-P. 


t  (E) 


where  we  have  defined  an  average  collision  time  t (E)  which 
is  a  function  of  the  electronic  energy  E. 

Equilibrium  requires 


LATTICE 


FIELD 


We  define  a  current  density  J(E)  by 


0 


3(E) 


e* 


(2) 


(3) 


(4) 
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where  m  is  the  effective  mass  of  an  electron  in  a  conduction  band. 
The  energy  given  up  by  the  electric  field  may  be  obtained  in  a 
straightforward  manner  from  the  Maxwell  equations  as  J*P.  Using 
Equations  (2)- (4)  we  then  have  for  the  net  energy  per  unit  volume 
transferred  from  the  field  to  the  electron 

A($,E)  =  •  (5) 

The  net  rate  of  energy  transfer  to  the  lattice  may  be  calculated 
by  summing  over  all  lattice  modes,  the  probability  that  the  ^ 
electron  will  either  absorb  or  emit  a  phonon  of  wave  vector  q. 


B  (E)  =  l  flu 

q  4 


emission 

(k-q,k) 


pabsorption 
F (k+q,k) 


(6) 


The  total  Hamiltonian  of  the  system  is  taken  to  be 

H  =  H  ,  +H  ,  +H  . 

electron  phonon  e-ph 

with  He_pk  representing  the  electron-phonon  interaction, 
considered  small  relative  to  the  total  Hamiltonian. 


The  quantum  mechanical  probabilities  of  emission  and  absorption 
.re  then  given  by  the  usual  "Golden  Rule"  transition  rate  formula 


Pe(k-q,k)  =  ^— | <k-q | Hg_ph | k> | 2 (nq+l) 


(8) 


P3 (k+q,k)  =  rL!<k+q!He_ph|k>|2nq6{|i(^+q)2  -  ft«q  -  (9) 


The  matrix  element  <&±q|  Hg  jJ  £>  represents  the  probability 

amplitude  for  an  electron  initially  in  a  state  o^  wave  vector 
|k>  to  scatter  to  a  final  state  of  wave  vector  |k±q>.  6  is  the 

Dirac  delta  function  ensuring  conservation  of  energy  and  n  is 
a  phonon  density  function  given  by  the  quantum  statistical4 
distribution  for  Bose  particles 


"<5  '  e*VkT-l 


(10) 
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The  condition  for  the  existence  of  a  steady  state  is  then  given 
by 


**^  (E),  .  |8.  |  fwq<  |  <k-q  |  He.ph  I  k>  | 2  (nq+l )  I , 


(11) 


- 1  <  Jc+q  I  He-ph  |  Jc>  I  2  nq<5  2 }  . 


6 j  and  6 j  refer  to  the  delta  functions  of  Equation  (8)  and  (9) 
respectively.  To  complete  the  calculation  it  is  necessary  to 
specify  functional  forms  for  J(E)  and  for  the  interaction 
matrix  elements. 


The  left  side  of  Equation  (2)  represents  the  rate  of  change 
of  momentum  pz  due  to  electron-phonon  interactions.  The  total 
rate  of  electronic  momentum  change  may  be  obtained  by  summing 
the  probabilities  that  the  electron  will  either  absorb  or  emit 
any  lattice  phonon.  We  may  thus  rewrite  Equation  (2)  in  the 
form 


1 

7(E) 


l  — TT"  l<k'<JIVph|k>l2(Vm' 

<k,+k> 

+  k  l<K+<jl>1e-ph|k>,2V2  • 


(12) 


2.  THE  ELECTRON-PHONON  INTERACTION. 

a.  Polar  Molecules.  The  electron-phonon  interaction 
term  will  depend  on  the  polarization  properties  of  the  molecular 
lattice.  In  polar  molecules  the  longitudinal  optical  modes 
produce  long  range  dipole  fields  with  which  the  electron  may 
be  strongly  coupled,  and  interaction  with  acoustic  phonons 
is  neglected. 

The  energy  of  an  electron  in  a  field  is  given  by  the 
expression 

He-ph*e*<r)  •  (13> 
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\£here  (r)  is  the  scalar  potential  function. 
D  is  related  to  the  polarization  field  ?  by 


The  divergence 


6  =  £+4n$  . 


The  vanishing  of  the  divergence  then  gives  the  polarization 
produced  by  the  optical  phonons  as 

?  -  w  • 

The  polarization  and  the  potential  are  expanded  in  fourier 
components 

I  vv*Svl5-S>  ' 


_  .  -+  -V  .  -+  + 

.  /v.v  _  1  r  .  _iq*r , ,  +  -iq»r  . 

-Mr)  -  jj*  l  ‘f’qe  H  +*q  e 


(14) 


(15) 


(16) 


(17) 


C  is  a  constant  which  will  be  determined  and  the  transformed 
coordinates  a  ,a+  correspond  to  the  phonon  annihilation  and 

creation  operators  in  second  quantized  representation. 
Substituting  (14)  and  (15)  into  (13)  we  obtain 


-4nic  _  .+  4nic  _+ 

♦q - q —  aq '  q - q~  aq 


(18) 


Substituting  (15)  and  (16)  into  (11)  gives 


H 


e-ph  /v 


.  -k.  I  42iS£  {a  e^^-a  +  e'1?'5)  • 

/\T^  L  a  K  a  a 


(19) 
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The  constant  C  is  calculated  in  Appendix  A  and  is  shown  to  be 


1  ,  hu> ,  *s  ,  1  1  , 

-  nr  It* j  “  t  j 


C  =  Iff  {3T 


e_  e 


Equation  (19)  gives  the  Frohlich  interaction  term  for  polar 
crystals . 

Using  (17)  and  the  selection  rules  for  harmonic  oscillater 
transitions  (cf.  Appendix  A)  we  have  for  the  absorption  and 
emission  transition  rates 


oa  -  2n 

p  -  fT 


e2  hu  1  1  ng<S  1 

“ LC"  $ 


pe  2H  e^_  fioi 
h  V  2 


ri_.ii 
Le»  eJ 


(V1) 


b.  Non-Polar  Molecules.  F.  Seitz  has  shown  that 
He-ph  arise  largely  from  the  electron-acoustic  phonon 

coupling  in  non-polar  molecules.  Let  V(r-x  ^)  equal  the 

potential  energy  of  an  electron  at  position  r  due  to  the  yth 
ion  in  the  th  unit  cell.  denotes  the  lattice 

vector  to  the  origin  of  the  unit  cell. 


X  .  =  £  +  u  .  , 

M  ,  l  V  r l 


where  u  represents  the  displacement  from  the  equilibrum 
u  i  * 

position  .  To  first  order  the  total  potential  at  r  may  then 
be  given  by 


*Frohlich,  H. ,  "Electrons  in  Metals,"  Adv.  in  Phys.  3,325  (1954) 
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V (r)  =  l  V (r-X  )  =  l  V(r-R  )  +  u  -7V(r-R  )  .  (24) 

1,V  i,v  l'u  W/1  w 

We  see  that  the  lattice  vibrations  induce  an  energy  shift  from 
the  value  perceived  by  an  electron  in  the  undistorted  lattice. 

This  energy  shift 


AV  =  l  u*  VV  (r-R  ) 


(25) 


is  considered  to  be  the^quantum  mechanical  perturbation  operator, 
where  the  displacement  u  &  is  a  phonon  operate*  and  vv  operates 

on  the  electronic  states  only. 

The  phonon  operator  u  is  expanded  in  normal  coordinates 
giving  w,e 

AV  =  T  ^V(r-R  .)  y  (—5 — )  *5e  ‘  Rw  £  (a+  -  a  )  ;  (26) 

nL  p,£  L  2MN(ua  '  e  a  a 

l  y  a  a 

a  represents  the  sum  over  both  wave-vector  and  polarization  and  M 
is  the  ionic  mass.  When  (26)  is  substituted  into  the  matrix 
element  for  electron-phonon  coupling  one  obtains  for  the 
transition  probabilities 


2n 

h 


I2hq  (n  +1) 5 1 


p-  =  ^. _ a_ 

y  *  2MsVn 


(27) 


..a  _  2n  l2«<3"q'2 
P  h  2MsVN 


(28) 


In  deriving  (27)  ,  (28)  we  used  the  acoustic-phonon  dispersion 
relation  (Dg  =  sq  with  s  the  sound  velocity.  Use  of**  Green's 

Theorem  allows  conversion  of  the  electron  state  matrix  element 
involving  vv  into  an  elect"';:-  overlap  integral  of  the  form 


I  =  D  /  '  d3r  . 


(29) 


TCI 


Ziman,  J.  M. ,  Electrons  and  Phonons,  o  c^rd  V.  Pres.,  (1960) 
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15 

Estimates  of  I  in  semiconductors  give  values  of  the  order  of 
one  electron  volt. 

3.  RELAXATION  TIMES  ^  The  term  (k^-k)/k  appearing  in  the 

expressior.  for  the  relaxation  time  (Equation  (12) )  may  be  obtained 
from  energy  and  momentum  conservation.  Conservation  of  energy 
gives 

cos(k,q>  «  ±  5jj|  ,  (30) 

where  the  minus  or  plus  sign  refers  to  emission  or  absorption, 
respectively.  Averaging  over  the  azimuth*’  angle  gives 


k  o 

z  _  -q*  -  nth oi 

kz  2k2  k2 


(31) 


We  assume  a  dense  distribution  of  modes  q  and  pass  from  a  sum 
to  an  integral  with  the  standard  approximation 


l 

q 


v 

(2nh) 3 


/  d3p 


(2nfi) 


n  /  p2dpd(cos  e) 


(32) 


\£here  £  =  hq  is  the  phonon  momentum  and  0  is  the  angle  between 
k  and  q,  the  initial  electron  momentum  and  phonon  wave  vector. 
The  integrations  over  the  probability  functions,  either  of 
the  form  of  (21-22)  or  (27-28),  involves  an  integral  over  a 

delta  function,  6  (^(kiq)  -  ■  -  ±  hiD^)  ,  as  defined  in 

(8-9).  The  integral  over  cos  9  can  then  be  written 


15Ershov,  A.  P. ,  "Ionization  During  Detonation  of  Solid  Explosives," 
Fixiba  Gor.  Vary,  11,  938  (1974). 
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/  a  (cos  6)«(&ia!  ±  i  tt£a  COS  0} 

*  /  m  rr  tn 


2m  " "q  m 

ilto  m 


m  . 2c2— 

_  f  d (cos  0)6  (cos  0  +  } 

“  J  "h  kq  2irihzhq  / 


(33) 


h2kq/m 


where  we  have  used  the  well  known  property  of  the  delta  function 


6  {  f (x) }  =  l 


6  (x  - 

1  / 


i  |f’  (x.) | 


(34) 


where  xi  are  the  roots  of  f (x) . 

The  integration  over  cos  0  then  gives 


/  d (cos0 )  6 (f  (0) )  = 


m 


h2kq 


(35) 


Equation  (12)  may  now  be  written 


Vm 


(/  q3G_(l  +  >  (n  +  l)dq 


tW  2(2u)2fi2k3 


■ftq2 


?r>o> 


+  J  *  V1  "  )nqdq  }  ' 


(36) 


where  is  just  the  probability  function  without  the  delta 
function  factor. 


For  electron  energies  E>>fiw  the  factors  l+2mu>/hq2  will 

4 

both  be  approximately  equal  to  one.  The  integral  may  then  be 
approximated  by 


4 


i 
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Ti&rmr&  I  '33°q(2nq  +  l)dq  . 


(37) 


For  long  wave  length  phonons  the  integrand  goes  to  zero.  We 
may  therefore  set  the  lower  limit  of  integration  equal  to  zero 
while  the  upper  limit  may  be  approximated  by 


g 

amax 


2k  . 


(38) 


For  polar  molecules  it  is  found  that  the  optical  phonons  are 
relatively  dispersionless.  The  factor  2nq+l  may  then  be  removed 
from  the  integral.  Using  (21-22)  in  (36)  we  have 


where 


1 

TTe] 


a  + 


2 

•haj/kT 

6 


a 


e2  /mai  (_i  _  i 

4tt/2  e°°  e' 


(39) 


(40) 


For  non  polar  molecules  we  use  Equation  (27-28)  and  approximate 


n  =  1  s  kT 

q  £V.sq/kT_  ^  -ftsq 


(41) 


The  relaxation  integral  then  takes  the  form 


,  /2 (m) 3/2I2k„T/E 

1  _  _ £  -  - 

TTeT  ~  2iTS2NMh4 


(42) 


where  it  has  been  assumed  that  the  main  acoustical  interaction  is 
with  the  long  wavelength  phonons. 
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4.  RATE  OF  ENERGY  TRANSFER  TO  THE  LATTICE.  The  energy 
transfer  rate  to  the  lattice  may  now  be  written  down  as 


B(E)  =  (2,)*fck  / 


(43) 


In  polar  crystals  the  integral  diverges  as  q->0.qm^n  is  then 

obtained  by  conservation  of  energy  (equation  along  with  the 
approximations ) 

k  •  k'  =  2k 


(44) 


k  -  k'  *  q  . 

^mm 


(45) 


This  gives  as  the  lower  limit  in  (40) 


moi 

3min  -fik 


(46) 


Performing  the  integration  we  then  have  the  following  expressions 
for  B(E).  In  polar  crystals 


B(E>  S'  . 


(47) 


In  non-polar  crystals,  taking  the  lower  limit  of  the  integral 
equal  to  zero  we  have 

D(E%  _  I2 (2m)3/2mE3/2 

B(E) - ?F?Fnm -  • 
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5.  CALCULATION  OF  CRITICAL  FIELDS .  The  critical  field  PB 
for  breakdown  is  now  calculated  from  the  energy  balance 
Equation  (11)  with  a  particular  choice  of  electron  energy  E. 

If  E  is  chosen  equal  to  the  ionization  energy  Ej  we  obtain  the 
following  expressions  for  the  critical  field  dependence  on  energy 
and  temperature.  For  polar  modes: 


3/2 


iTTur 


meio 


(—  -  -) 
£  CO  £ 


/El 


4Et 

{1°9  ET 


}  U  + 


3fia)/kT_1J 


(48) 


For  acoustic  modes: 


2f 

zWs  M  1 


(49) 


In  Figures  1  and  2,  we  have  graphed  the  behavior  of  the 
critical  field  as  a  function  of  temperature  for  a  reasonable 
choice  of  physical  parameters. 

To  illustrate  the  use  of  Equation  (48)  for  a  typical  polar 
molecule,  we  will  estimate  the  breakdown  voltage  for  TNT.  In 
Figure  1,  we  have  chosen  E^  as  6e.v.l5  and  have  taken  for  a> 

a  typical  frequency  of  10^4.  The  other  parameters  chosen  were 

T=300°K,  m=m  (electron)  and  (l/e~-l/e)  =~— where  eQ  is  the 

ze0 

free-space  permittivity.  We  then  obtain  for  the  critical  breakdown 
field,  Fc-1()6  V/cm.  This  value  should  be  accurate  to  within  an 
order  of  magnitude  given  the  uncertainty  in  the  physical  parameters. 

6.  TIME  LAG  BEFORE  BREAKDOWN.  We  shall  make  a  rough 
estimate  of  the  time  till  breakdown  assuming  a  critical  field 
Fb  =  10 6  v/cm  has  been  applied. 

Let  the  electron  be  accelerated  from  rest  and  assume  that 
it  ionizes  another  electron  immediately  upon  attainment  of  the 
critical  energy 


^mVg  =  I 


(50) 
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If  we  take  I  =  6e.v,  characteristic  of  TNT  then 

VB  =  =  1.4  x  106m/s  .  (51) 

-14 

With  a  relaxation  time  t  =  10  sec,  an  electron  will  achieve 
a  velocity 

eF 

Vn  =  — -  t  =  1.7  x  105m/s  (52) 

before  collision,  which  is  not  enough  to  achieve  ionization.  We 
must  then  ask  how  many  collisions,  on  the  average,  will  occur 
before  an  electron  will  have  sufficient  time  tg  between  collisions 
to  achieve  the  critical  velocity  Vg.  We  can  estimate  this  by 
writing  down  a  probability  function  P(t)  representing  the 
probability  that  a  time  t  has  passed  with  no  collisions.  We 
shall  take  a  model  P(t)  as 


P(t)  =  e~tA  ,  (53) 

which  has  the  proper  boundary  behavior 

P(o)  =1  and  P(“)  =  0  ,  (54) 

and  which  also  falls  off  rapidly  for  t>T.  We  may  also  formally 
derive  Equation  (53)  by  noting  that  P(t+dt)  represents  the 
probability  that  no  collision  will  occur  in  a  time  t  or  in  a 
following  time  interval  dt.  We  thus  have 

(55) 

boundary  condition 

(56) 

formally  obtain 


P  (t+dt)  =  P(t)  P (dt) 


expanding  P(dt)  about  t  =  0  and  applying  the 
we  have 


P (t+dt)  =  P(t) 

If  we  make  the  identity  x  ^ 
the  equation 


f1+ 


dP 

It 


-dP 

It 


t=0 


we  may 
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dP  _  -P{t) 

HIT  t  .  (57) 

Integrating  (57)  gives  us  the  probability  function  of  Equation  (53) . 

The  time  needed  for  a  single  electron  to  achieve  breakdown 
is  given  by 

tQ  =  mVg/eFg  5  9  x  10  ^  sec  (58) 

with  the  values  of  Vg  and  Fg  given  above.  We  thus  have 

P(tfi)  =  e~9  5  1/8000  .  (59) 

It  will  take  on  the  average  8000  collisions  to  produce  a  single 
ionization  event.  The  time  for  each  ionization  event  is  thus 
given  by 

tj  *  8000t  =  8  x  10-11  sec  .  (60) 

To  estimate  now  the  total  number  of  ionizing  collisions  n  which 
will  characterize  dielectric  breakdown  of  the  material  we  equate 
the  energy  supplied  by  the  field  to  the  electrons  with  the  energy 
required  to  raise  the  lattice  to  some  critical  temperature  T0. 

After  n  collisions  we  have  2n  free  electrons  so  we  set 


2nI  =  cvpV(To-T)  .  (61) 

V  represents  the  volume  of  the  cone  to  which  the  2n  electrons 
are  assumed  confined.  The  breakdown  time  Tg  will  be  given  by 


Tfi  =  ntj.  .  (62) 

The  volume  of  the  cone  V  is  given  by 

V  =  1/3JI  (vig) 3tan2a  (63) 

where  v  is  the  average  velocity  in  the  field  direction  given  by 
the  condition 
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v  - 


^5  V, 


B 


(64) 


and  a  is  the  half  angle  of  the  cone.  Tan  o  may  be  evaluated  by 
assuming  that  motion  perpendicular  to  the  applied  field  is  governed 


only  by  thermal  diffusion, 
by 


The  thermal  velocity  V.  is  determined 


„t .  (*?)' 


(65) 


After  a  time  tb  corresponding  to  n  collisions  the  electrons  may 
diffuse  out  to  a  distance 


i 


n 


*5 


2kT  ^ 
m  B 


(66) 


where  the  n  term  comes  from  the  drunkard's  walk  solution.  We 
then  have 


tan  a  =  (2knT) 


(67) 


Substituting  (62-67)  into  (61)  and  taking  the  log  of  both  sides, 
we  obtain  a  transcendental  equation  for  n 


4  log  n 
log 


cyp  ATirkgT 
31m 


log  2 


(68) 


We  take  cv  =  .2  cal/gm,  p  =  1.6  gm/cm3,  AT  =  100°K  with  the  other 
constants  as  before.  We  have  solved  Equation  (8)  graphically  («ee 
Figure  3)  and  obtain  for  n 


n  *  52 


(69) 


We  then  have 

■tg  =  52tB  =  4  x  10  ®  sec  .  (70) 

Thus  the  time  delay  until  breakdown  after  the  sudden  application 
of  a  critical  field  is  of  the  order  of  nanoseconds. 
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7.  CHARACTERISTIC  RESULTS.  The  prediction  of  actual 
breakdown "fields  ?rom  formulas  such  as  (48)  and  (49)  requires 
knowledge  of  various  material  parameters  such  as  the  effective 
mass,  dielectric  constants,  ionization  potentials,  etc. 

Unfortunately  many  or  most  of  these  have  not  been  measured  for 
explosive  materials.  Lack  of  an  adequate  experimental  data 
base  precludes  accurate  predictions  of  breakdown  strength  in 
explosives.  Some  general  characteristics  of  breakdown  may, 
however,  be  given.  These  include: 

1.  It  occurs  very  rapidly.  Experimental  breakdowns 

—  8 

are  observed  to  occur  in  less  than  10  sec,  with  some  instances 

,.-6 

as  rong  as  10  awx. 

2.  It  is  the  breakdown  of  an  equilibrium  condition  and 
therefore  occurs  at  an  essentially  constant  temperature.  Breakdown 
may  thus  occur  at  low  or  room  temperatures. 

3.  In  the  type  of  breakdown  we  have  considered,  dielectric 
strength  will  be  independent  of  sample  geometry. 

Other  modes  of  breakdown  also  exist.  In  unpublished  work 
recently  completed  at  NSWC  J.  Forbes  measured  breakdown  fields 
in  a  number  of  explosives.  In  the  experimental  setup  employed, 
the  explosive  sample  was  sandwiched  between  two  electrodes. 

Breakdown  theory  as  developed  in  the  previous  section  would  not 
be  applicable  as  we  would  expect  substantial  carrier  injection 
from  the  electrodes^® .  This  would  lead  to  lower  measured  critical 
fields  than  might  be  expected  in  an  "electrodeless"  geometry. 

8.  CRITICAL  ELECTRON  DENSITY.  When  the  density  of  electrons 
is  large  enough,  the  single  electron-phonon  interaction  picture  is 
no  longer  realistic  and  the  effect  of  electron-electron  collisions 
must  be  included  in  any  calculation. 

The  critical  density  for  the  validity  of  the  one-electron 
approximation  is  obtained  by  considering  the  net  rate  of  energy 
transfer  to  the  lattice.  As  the  rates  of  absorption  and  emission 
of  a  lattice  phonon  are  given  by  nq  and  (1+nq)  respectively, 

it  will  take,  on  the  average  2nq+l  collisions  to  transfer  a  single 

quantum  of  energy  hw.  The  transfer  rate  may  then  be  approximated 
by 


dU  _  1  ha' 

dtLATTICE  (2nq  +  H  tTeT 

^6Fowler,  R. ,  and  Nordheim,  L. ,  "Electron  Emission  in  Intense 
Electric  Fields,"  Proc.  Roy.  Soc.,  A119 ,  173  (1928). 


(71) 
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The  rate  of  energy  loss  due  to  collision  with  other  electrons  is 
given  by^7  Bohm  and  Pines  as 


dU 

d “ELECTRON 


47tne‘4 

/2mE 


(72) 


The  critical  electron  density  for  a  consideration  of  electron- 
electron  effects  is  obtained  when  (71)  and  (72)  are  of  equal 
magnitude.  This  gives  the  condition 


n 


CRITICAL 


rtfiE 

4ire4 


hw 


(2nq  +  1)  tTE) 


(73) 


In  polar  materials,  substitution  of  some  representative  values 
gives  nCRITICAL  =  1017(cm)  3  typical  values  for  non-polar  crystals 


gives  nCRITICAL  =  lO14'15 (cm) -3. 


18 

Published  values  for  conductivities  in  explosive  materials 
give  values  of  a  in  the  range  of  10~^°  -  10“8  (ft-M)"l.  Associated 
number  densities  are  obtained  from 


n  = 


ma 
e2  r 


(74) 


Taking  the  effective  mass  equal  to  the  electron  mass  and  using 

and  use  of  the 
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a  value  of  x  of  10  we  obtain  n<n, 


CRITICAL 

single  electron  approximation  is  therfore  justified. 

Recently  G.  Hammond  (NSWC)  has  undertaken  some  conductivity 
measurements  of  explosives  in  the  microwave  regime.  Preliminary 
results  indicate  substantially  enhanced  conductivities  at  these 
frequencies.  We  shall  therefore  briefly  outline  the  modifications 
to  the  single  electron  treatment  given  above  which  enable  treatment 
of  the  many  body  electron  distribution. 


17Pines,  D. ,  "A  Collective  Description  of  Electron  Interactions," 
Phys.  Rev.  92,626  (1953). 

18Federoff  and  Sheffield,  Ed.  Encyclopedia  of  Explosives  D1221 
(1972) . 
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9.  EFFECTS  OF  THE  ELECTRON  DISTRIBUTION.  When  the  many 
electron  system  is  describee!  by  a  distribution  function  f  defined 

2  3  3 

in  such  a  way  that  fd  qd  p  describes  the  number  of  electrons 

X2Trh)3  n 

m  volume  element  dJa  with  wave  vector  k=e,  the  current  density 
may  no  longer  be  given  by  Equation  (4)  but  is  now  defined  by 

J  =  |  /  vfd3k  .  (75) 


In  general,  only  the  asymmetrical  component  of  f  will  give  rise  to 
a  current.  Equations  such  as  (6)  and  (12)  representing  probabil¬ 
ities  for  electronic  transitions  must  be  augmented  by  factors  such 
as  f (k) ( 1— f (k' ) )  which  represent  the  probability  f(k)  that  the 
initial  state  | k>  is  occupied  and  the  final  state  |k’>  is  unoccupied. 
Following  this  prescription,  rhe  current  density  J  and  the  energy 
transfer  rate  B  may  be  computed  and  a  breakdown  field  Fc  derived 
as  in  section  (5)  which,  however,  accounts  for  the  effects  of 
the  electron  distribution.  It  has  been  implicitly  assumed  here 
that  the  electrons  and  phonons  are  thermally  equilibriated. 

In  the  high  density  case,  the  assumption  of  an  electron-phonon 
system  in  thermal  equilibrium  may  no  longer  be  valid.  With  a 
sufficient  density  of  conduction  electrons  the  electrons  may 
exchange  energy  among  themselves  at  a  much  greater  rate  than  with 
the  lattice.  The  phonons  and  electrons  must  then  be  treated  as  a 
two-fluid  system  with  the  phonons  at  temperature  T  and  electrons 
at  temperature  Te.  The  applied  field  imparts  momentum  to  the 
electrons  shifting  them  away  from  the  0-field  equilibrium  momentum 
to  some  new  p.  This  particular  problem  has  been  treated  extensively 
by  Frohlich  and  Paranjape3  (1956)  and  we  shall  outline  their 
calculation  here.  The  equilibrium  condition  is  given  by  the 
vanishing  of  the  Boltzmann  equation 


^Lttice  +  (“Leld  + 


=  0 


(76) 


FIELD 


'ELECTRON 
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thatrSeehaveni°  collisions  conserve  both  energy  and  momentum  so 
I  p(ff)  =0  and  J  El  III)  .  Q 

P  V3t'e  p  2m  V 3 1  /0  "  0  '  (77) 

Equation  (55)  then  gives  the  two  conditions 

l  «(&  *  (Mb  ■ . 


5*88  •(«)•-.. 


P  . .  L  *  “  w  F 

The  distribution  function  is  assumed  to  be  of  the  form 
f(p)  =  ae~(P  -  P0)2/2mkTe 


(78) 

(79) 


(80) 


The  fcw°  Equations  (78)  and  (79)  are  then  4.  . 

to  determine  the  two  unknowns  d  and  S  sufficient,  m  principle, 
computation  of  '  a  mean  coStioS  Snslty  Calculation  Po  enables 


nep 

J  = 

m 


(81) 


The  rate  of  energy  transfer  to  the  lattice  will  be  given  by 

B<V*>  -/£(!!)  ' 


(82) 


where  the  summation  in  Equation  (79)  h  =  ~ 

intearal  "  (79)  has  been  converted  to  an 

3t ^LATTICE  is  ofc)tained  formally  by  consjderinq  the 
probabilities  of  transitions  into  or  out  of  stales  of  llZZ,. 
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=  l  f(p)Pa(p,p  +  q)  -  f  (p  +  q)PS(p  +  q,p) 
L  q 

+  f(p)Pe(p,p  -  q)  -  f(p  +  q)Pa(p  -  q,p)  . 


(83) 


It  is  found  that  Equation  (82)  has  a  maximum  at  some  Te.  Setting 
J‘F  equal  to  B(T,TeMAX)  then  defines  a  critical  breakdown  field. 

10.  APPLICATION  TO  EXPLOSIVES.  The  application  of  electrical 
breakdown  theory  to an  explosive  medium  follows  closely  breakdown 
theroy  in  inerts.  Since  electrical  breakdown  is  perceived  as 
the  destruction  of  steady  state-equilibrium  condition,  the  most 
important  consideration  is  the  nature  of  the  electron-phonon 
coupling.  Important  parameters  are  the  temperature  T  the  phonon 
spectrum  and  the  various  optical  constants.  These  considerations 
apply  equally  well  to  explosives  in  metastable  equilibrium  with 
their  environment  as  of  inerts.  It  is  only  in  a  consideration  of 
thermal  breakdown  that  the  exothermic  kinetics  of  explosive 
materials  must  be  accounted  for. 


III.  THERMAL  BREAKDOWN 

1.  THERMAL  BREAKDOWN  IN  INERTS.  In  thermal  breakdown  it  is 
assumed  that  the  properties  of  the  lattice  may  be  fully  described 
by  the  macroscopic  temperature  dependent  conductivity.  As  an 
electric  field  is  applied,  energy  will  be  deposited  through  Joule 
heating.  This  deposited  energy  will  raise  the  temperature  of 
the  lattice  which  results  in  an  increase  in  the  conductivity 
function.  The  increased  conductivity  gives  rise  to  an  increased 
current  which  implies  a  higher  rate  of  joule  heating  and  so  on. 

If  the  lattice  is  unable  to  dissipate  this  energy  at  a  rate  great 
enough  to  overcome  the  joule  heating,  temperatures  will  rise 
uncontrollably  and  thermal  breakdown  will  be  said  to  occur. 

The  critical  field  is  then  defined  by  the  greatest  field  for  which 
a  steady  state  solution  of  the  heat  transfer  equation  (Equation 
(84))  is  possible.  Formally,  the  theory  of  thermal  breakdown 
in  inerts  consists  of  finding  solutions  to  the  equation 

V(KVT)  +  a  (F,T) F2  (t)  .  (84) 
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cv  is  the  specific  heat  per  unit  mass,  p  is  the  density,  k  is 
the  thermal  conductivity .  a represents  the  rate  of  energy 
deposited  through  Joule  heating.  Equation  (84)  says  that  the 
input  energy  is  partly  stored  as  internal  energy  tending  to 
increase  the  temperature  and  partly  conducted  away. 

Many  dielectrics  have  been  found  to  follow  a  Poole-Frenkel 
conductivity-*-9  law  given  by 


<J(F,T)  =  aQ(F)e 


-*/kt 


(85) 


<f>  represents  a  work  function  for  electron  ionization  from  a 
valence  band  or  a  trap  to  the  conduction  band,  a  may  in  general 
be  a  function  of  the  electric  field  F. 

2.  THERMAL  BREAKDOWN  IN  REACTIVE  MEDIA.  With  the 
identification  of  the  Poole-Frenkel  conductivity  function  the 
theory  of  thermal  breakdown  in  inerts  is  seen  to  be  mathematically 
equivalent  to  the  theory  of  thermal  explosions  with  the  joule 
heating  term  a F2  taking  the  place  of  the  usual  Arrhenius  kinetic 
term  in  thermal  explosion  theory.  If  we  now  wish  to  consider 
thermal  dielectric  breakdown  in  explosive  materials,  Equation  (84) 
must  be  augmented  by  an  appropriate  heat  source  term  to  account 
for  the  chemical  reactivity  of  the  medium.  The  appropriate 
equation  for  thermal  breakdown  in  explosives  is  now  given  by 


c  o|~r  =  F2c  +  pQZe“Ea/kT  +  V(ieVT) 

V  a  u  O 


(86) 


Equation  (86)  can  equally  well  be  though  of  as  representing  the 
influence  of  an  applied  electric  field  on  thermal  explosion 
solutions.  An  applied  electric  field  will  influence  thermal 
explosion  parameters,  e.g.  induction  times  to  explosion,  only 
insofar  as  the  termC^F2exp  (— 4>kT)  is  not  negligible  as  compared 
to  the  Arrhenius  kinetic  heat  generation  term.  Alternatively, 
the  concept  of  thermal  breakdown  of  an  explosive  remains 
meaningful  only  so  long  as  the  J*F  term  is  not  overwhelmed  by 
the  Arrhenius  term. 

_________________________ 

Frenkel,  J. ,  "On  Pre-Breakdown  Phenomena  in  Insulators  and 

Electronic  Semiconductors,"  Phys.  Rev.  54,  647  (1938). 
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As  is  well  known,  analytic  solutions  to  Equation  (86)  may 
not  be  written  down  in  general  because^  of  the  highly  non¬ 
linear  source  terms.  A  variety  of  approximate^'^  solution 
schemes  for  various  geometries  have  been  given  in  the  literature. 
It  is  also  possible  to  solve  (86)  in  various  limiting  cases. 

These  include  the  steady  state  regime  obtained  by  setting^/ 24 
3T 


at 


=  0  and  the  "impulse  thermal  breakdown"  obtained  by  setting 


V‘(kVT)  =  0.  We  have  chosen  to  solve  Equation  (86)  numerically 
retaining  all  terms  except  the  surface  loss  terms,  assumed  to 
be  small  for  the  boundary  conditions  chosen. 


3.  CRANK-NICHOLSON  METHOD.  Equation  (86)  is  an  example  of 
a  parabolic  partial  differential  equation  with  a  non-linear  source 
term.  The  Crank-Nicholson^S  scheme  provides  a  numerical  algorithm 
for  evaluating  (86)  which  is  correct  to  second  order  in  time  and 
stable  for  any  values  of  increments  Ax  and  At.  Figure  4  represents 
the  space-time  grid  for  the  numerical  calculation.  At  time  t=0 
the  solution  is  given  at  all  spatial  points  i.  The  crank-Nicholson 
scheme  involves  writing  finite  difference  approximations  for  the 
differential  operators  at  points  x^,tn+^  where  tn+^  represents 

points  halfway  between  the  known  level  at  n  and  the  unknown  time 
level  n+1. 

Equation  (86)  is  solved  in  a  slab  geometry  illustrated  in 
Figure  5  with  boundaries  at  x=0  and  x=xf.  The  divergence  term 
then  takes  on  the  form  of  a  second  derivative,  the  thermal  conduc¬ 
tivity  <  assumed  constant.  The  finite  difference  approximations 
to  the  derivatives  are  given  by 


3T 

at 


i,n+J$ 


T  -  T 

i,n+l  i,n 

At 


(87) 
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FIGURE  4  GRID  FOR  FINITE  DIFFERENCE  COMPUTATION  IN  A  CRANK-NICHOLSON  SCHEME 


i 
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fih\ 

v«!A, 


n+*5 


T i+l,n  ~  2Ti,n  *  Ti-l,n 
(Ax)2 

+  T. 


(88) 


+  Ti+l,n+l  2Ti,n+l  n  ii-l,n+l 
(ax)  2 

The  source  term  is  an  explicit  function  of  temperature  and  is 
estimated  by 

■  *  [S(Ti,„>  +  S(Ti,n+l>]  • 

S  (T^  n+1)  will  be  given  by  the  approximation 

s<Ti,n+l>  "  S(T*>  +  (ft).  _  <Ti,n+l  -  T*>  ' 


i,n 


where 


T*  "  2Ti,n  -  ~i,n-l  ' 


(89) 


(90) 


(91) 


The  algorithm  for  the  temperatures  at  the  unknown  timelevel  n+1 
in  terms  of  known  values  at  time  level  n  is  then  obtained  by 

substituting  (87-91)  into  (86) .  We  then  have 

/ 

m  +  Lo  _  2  (AX)  2  as  (Ax)  2  I 

i-l, n+1  \  A  KAt  3 T*  2k  j  i,n+l  i+l,n+l 

=  -T.  .  +  <2  -  ---rr- /  T.  -  -  T.  .  -  h  <  S (T .  )  +  S(T*) 

i-l, n  |  KAt  |  i,n  i+l, n  |  i,n 

-  T(lf)T.J  . 

The  boundary  conditions  are  obtained  letting  the  sample  sit  in 
a  heat  bath  at  temperature  T0.  We  then  have  Tx=TxFINAjj=To  ^or 

all  n.  Equation  (92)  defines  a  set  of  linear  equations  in  the 
unknown  timestep  n+1. 


(92) 
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Solutions  for  all  the  at  a  given  timestep  are  obtained 
simultaneously  by  inverting  the  coefficient  matrix  defined  by 
(92).  In  Appendix  B,  we  display  the  computer  program  written 
for  this  problem.  As  is  evident  from  Equation  (92)  the  coeffi¬ 
cient  matrix  is  in  tridiagonal  form.  All  programs  were  run  on 
the  NSWC/WOL  CDC  6500  computer.  In  Appendix  C  we  explicitly 
exhibit  the  tridiagonal  matrix  and  coefficients. 

4.  DISCUSSION  OF  RESULTS  OF  NUMERICAL  INTEGRATION.  In  Figure 
6,  we  show  the  solutions  generated  for  a  typical  set  of  explosive 

parameters.  We  have  taken  p=1.84  gm/cm1,  k  =  10  1  Q=1000  cal/gm, 

Z=10“14  (sec)"1,  C  =.27  ,  and  E  =2  e.v.(l  e.v.  3  20,000  cal/gm). 

v  _  gm°K  .  a 

We  have  taken  ao=10-/  (0-m)-1  in  line  with  Picatinney  data  and 
have  taken  the  electric  field  to  be  a  constant  during  each 
calculation.  The  constant  chosen  ranged  from  0  through  10°  V/cm. 

The  work  function  was  chosen  at  3  e.v.  The  general  results 
are  typical  of  thermal  explosion  calculations.  The  temperature 
in  the  sample  rises  relatively  quickly  to  the  bath  temperature  T0 
or  slightly  above.  If  T0  is  chosen  below  some  critical  temperature, 
computed  here  as  540°K,  the  temperature  stabilizes  at  T0.  When  T0 
is  greater  than  the  critical  temperature,  the  temperature  rises 
initially  to  T0,  remains  there  for  some  characteristic  induction 
time  and  then  increases  very  steeply  -  indica-ing  thermal  explosion. 
There  is  no  indication  of  electric  field  dependence  in  Figure  6 
since  it  was  found  that  for  the  above  choice  of  parameters  the 
applied  electric  field  has  absolutely  no  effect  on  the  thermal 
explosion  characteristics.  It  is  then  clear  that  the  concept  of 
thermal  breakdown  of  an  explosive  loses  validity,  being  overwhelmed 
by  the  chemical  production  of  heat  and  subsequent  thermal  explosion 
at  a  time  well  before  joule  heating  alone  could  raise  the  temperature 
significantly. 

5.  INFLUENCE  OF  THE  CONDUCTIVITY  WORK  FUNCTION.  The  critical 
parameter  for  possible  electric  field  influence  is  the  work  function 
4> ..  With  <t>  chosen  as  3  o.v.  the  joule  heating  is  relatively 
inefficient  even  for  much  larger  values  of  ctq  corresponding 

to  Hammond's  preliminary  results.  Since  there  are  no  available 
experimental  measurements,  it  seems  most  reasonable  to  select  <J> 
to  be  somewhat  larger  than  the  activation  energy  Ea.  The  reason 
for  this  is  that  represents  the  energy  to  ionize  an  electron 
into  a  conduction  band  while  Ea  only  has  to  supply  enough  energy 
to  excite  the  molecule  to  the  "activated  complex"  energy  state. 

It  is,  however,  conceivable  that  under  some  ci reruns tances  $  might 
be  less  than  E  .  This  might  arise  if  the  conduction  mechanism 
did  not  arise  from  a  valence  to  conduction  band  transition. 
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PARAMETERS:  p  =  1.84  gm/em3 ,  z  -  1014  SEC  ,  Q  *  |1kcal,  E#  -  2#.v,  C„  -  .Seel/0!*, 

&  =>  1.25  x  106V/cm,  0  =  10  *6|  [V/cml*56 .  oQ  =>  10  '7  (J^-m)'1 ,  £  -  1cm. 

FIGURE  7  TEMPERATURE  EVOLUTION  FOR  Dl  FFERENT  BATH  TEMPERATURE  T„ 

D 
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One  possibility  which  suggests  itself  is  a  conduction  due  to 
the  ionization  of  trapped  impurities^ .  The  trapping  level  could 
lie  relatively  close  to  the  conduction  band  and  would  thus  require 
a  smaller  work  function.  In  Figures  7  and  8,  we  have  graphed 
solutions  indicating  the  possibility  of  radically  altering  the 
thermal  explosion  parameters  by  applying  a  strong,  step  function 
electric  field.  The  parameters  are  the  same  as  in  Figure  6. 

The  possible  tuning  of  explosion  parameters  through  controlled 
doping  of  the  explosive  sample  also  suggests  itself. 

G.  THERMAL  BREAKDOWN  CHARACTERISTICS.  The  characteristics 
of  thermal  breakdown  of  dielectrics  have  been  seen  to  be  indis¬ 
tinguishable  from  thermal  explosion  characteristics.  These  include 

a.  relatively  long  induction  times  until  explosion 

b.  dependence  on  geometry 

c.  high  temperature  phenomenon 

Unlike  electrical  breakdown  which  occurs  at  a  constant 
temperature,  thermal  breakdown  tracks  the  build  up  of  temperature 
during  the  breakdown  process. 

IV.  SUMMARY  AND  CONCLUSIONS 

We  have  analyzed  the  application  of  dielectric  breakdown 
theory  to  explosive  dielectrics.  In  a  lower  temperature  regime, 
because  the  explosive  is  in  metastable  equilibrium  with  its 
environment,  breakdown  theory  could  be  modeled  after  that  in 
inerts.  We  have  predicted  electric  breakdown  fields  on  the  order 
of  1  MV/cm. 

Ir.  a  higher  temperature  regime,  thermal  breakdown  theory 
had  to  be  modified  to  account  for  the  exothermic  kinetics  of 
a  reacting  medium.  It  was  concluded  that  the  distinction  between 
thermal  breakdown  and  thermal  explosion  disappears  for  low  enough 
values  of  the  conductivity  work  function.  At  still  smaller 
values  of  $  a  possible  radical  effect  of  applied  electric  field 
on  thermal  explosion  parameters  exists. 

A  resolution  of  these  questions  and  more  accurate  predictions 
of  breakdown  strengths  requires  a  substantial  expansion  of  our 
data  base  on  explosive  materials.  In  particular,  the  absorption 
spectra  of  explosive  molecules  through  the  infrared  is  of  fundamental 
importance  and  should  be  mapped  out.  The  various  material  constants 
such  as  the  effective  mass  or  the  high  and  low  frequency  dielectric 
constants  could  be  determined  e.g.  by  Raman  or  Brillouin  scattering 
experiments.  Knowledge  of  absorption  spectra  would  also  be 
prerequisite  to  an  experimental  program  designed  to  measure  molecular 

26 

Frohlich,  H.,  "On  the  Theory  of  Breakdown  in  Solid  Dielectronics, “ 
Proc.  Roy.  Soc.  A188,  521  (1949). 
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energy  elaxation  through  the  newly  developed  techniques  of 
picosecond  spectroscopy^.  More  measurements  of  the  conductivitie 
of  explosive  dielectrics  are  needed.  These  measurements  should 
be  made  at  various  temperatures  to  test  the  Poole-Frenkel  model 
and  determine  activation  energies  for  conduction.  Ionization 
energies  for  explosive  molecules  should  be  measured.  To  separate 
out  space  charge  and  electron  tunneling  effects,  measurements 
of  electrical  breakdown  should  be  made  with  electrodes  at  some 
distance  from  the  explosive  surface  in  an  evacuated  chamber. 
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APPENDIX  A 

CALCULATION  OF  THE  INTERACTION  CONSTANT 


To  determine  the  constant  C,  we  consider  the  interaction  energy 
of  two  point  charges.  Classically  this  energy  is  given  by 


E  = 


.Jr12c 

9 


(A-l) 


where  r.  represents  the  interparticle  distance.  The  dielectric 

function  e  represents  the  total  response  of  the  lattice  to  the 
perturbing  radiation.  This  dielectric  response  is  given  by 


e 


e  + 
00 


,  2 
(e<feJ  % 


(A- 2) 


represents  the  high  frequency  background  response  due  to  higher 
energy  electronic  resonances.  If  the  polarization  due  to  lattice 
displacement  (i.e.,  phonon  coupling)  is  ignored  the  interaction 
energy  would  be 


E 


(A- 3) 


The  net  energy  due  to  electron-phonon  interaction  can  therefore 
be  represented  by 


AE 


e 


4nr 

1/2 


1 

F 


(A- 4) 
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We  may  also  write  down  the  interaction  energy  shift  due  to  electron 
phonon  interactions  from  the  quantum  mechanical  perturbation  formula 


4E  =  I  l<n|H  h(r,)  +  H  (r2)|0>| 


He-ph  given  by  Equation  (19)  and  involves  a  sum  over  one- 

phonon  creation  and  annihilation  operators.  The  matrix  element 
is  given  explicitly  by 


(A-5 ) 


<n| H (r) | o>  = 


Ce4ni 


/V  q  q 


I<q|2iq-raa  -  e-iq'V|0> 


(A-6) 


Applying  the  selection  rules 


<n ' | a  | n>  =  /n+1  6  . 

1  q1  n+l,n 

<n' | a  +| n>=  /n  5 

1  q  1  n-l,n 


(A- 7) 


(A-8) 


we  have 


int  /v  <3 


(A- 9) 


The  denominator  of  (A-5),  E  -E  ,  represents  a  single  phonon 
quantum.  It  has  been  found  that  polar  optical  phonons  display 
little  dispersion.  The  denominator  may  then  be  removed  from  the 
sum  as  approximately  constant  and  equal  to  hw,  with  u  the  polar 
optical  frequency. 

Applying  (A-9)  in  (A-6)  and  neglecting  the  self  energy  term  we 
have  for  AE 


_  -2  (rile)  2C2  v  e1^-^ 

vHu  £  T2  • 


(A-10) 


A- 2 


NSWC  TR  79-331 


oooooo  nonnoi » 


NSWC  TR  79-331 


"'lOjuA*' 


c 


APPENDIX  B 

COMPUTER  PROGRAM  FOR  THERMAL  BALANCE  EQUATION 
UTILIZING  A  CRANK-NICHOLSON  SCHEME 


"IM  I  / .  '  /  3 


M  N  4  .  b  +  4  <>* 


PROGRAM  ^I<cl(IV;>jrt0UT3JTfTAcr«>  =  3UTD  JT) 


SET  IvITlAL  PARAMETERS 

XL  =  SLA3  THlCKVESSf  IMAX=  NJM-'ER 


L5VEL  \,T:  =  THE*<  vAL  SONOo'CTI  VJTYtRC=D£NSl  TY.|  =  M-  au?VCY  t 
ILAST  =  I  MAX-1 

1  n  =  i' a  x - ? 

M  T  =  5  j 
DT=.002L 

0X=XL/I max 

t:=i .£-? 
ro=i  ,e<i 

Z  =  1  .  E  1  4 
r-  =  l.E3 
OH  =  .  >  7 

X  x  r  1 

X;<  =  TC/(R0*SH) 

£1=2. 

E?  =  3. 

SI GMA=1 -1 

DELTA=.0i 

E  =  1  .  ?  5 E  •• 

INITIAL  C  0  V  )  I T  I  U  U  S  »  0  T  ,H  j  0  \  0  A  R  I F  $  F  '-LAI'  a  R  ~  AT  THT 

BATH  TEMPERATURE  UG»  INTERIOR  INITIALLY  AT  303 "<rL YIN 
U(I)  ALWAYS  RERE3ENTS  TrM=rRATJRES  AT  TnF  KNOW  TTHr'  Ki/'i  'i 
UP(I)  ARE  THE  UNKNOWN  TCMP?SATJ?r|  *{  t  ]  « f  L  r  G  r  ,  T  { Jf  Ljt<M  ‘ 
IS  AN  INTERPOLATION'  3ETWEM'N  A\5n+1.  U  J  5<I) 

U0  =  5  55 . 

U( IHAX>  =  JO 

,,  00  13  1=1,ILA3T 

UU)  =  30  3  . 
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BOUNDARY  CONDITIONS 
US  (I  '•AX  )  =U0 
LPt  IXAy>-jQ 
WRITE (S,S> 

b  c0P,MS  TU  wl  t/f  *  =AP  AMETER  SET*,/) 
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♦»*  ACTIVATION  f.NERSY  =  *,E3.?*/*  »  ORK  F  UNC  T  I  0  N  =  *  ,  E  9 . 2  ,  *  S13MA=**E9. 
♦♦*  TOTAL  LE3\7H<:'*.)s*9!r,-.l*/*  DE LX  =  * ,E9. 2  *  *  NO.  Oc  T IMCSTEPS  =  * ,  I 
OFLT-*,E9.2  ,/*  ELECTRIC  F I  EL  D  (  V/C  » >  =*  *  E  9.  1 ,  / //> 
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APPENDIX  C 

THE  TRI DIAGONAL  MATRIX  AND  THOMAS  ALGORITHM 


Equation  (71)  is  of  the  form 


aiui-i  +  Vi  *  Vi+1  -  di  < 

where  the  u^  represent  the  unknown  temperatures  at  the  spatial 
points  i  at  time  level  n+1.  The  coefficients  are  equal  to 
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(C-4) 


The  tridiagonal  matrix  equation  for  each  n  then  takes  the  form 
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